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Abstract 

We present a rigorous analysis of the phenomenon of decoherence for general 
N— level systems coupled to reservoirs. The latter are described by free massless 
bosonic fields. We apply our general results to the specific cases of the qubit and 
the quantum register. We compare our results with the explicitly solvable case of 
systems whose interaction with the environment does not allow for energy exchange 
(non-demolition, or energy conserving interactions). We suggest a new approach 
which applies to a wide variety of systems which are not explicitly solvable. 

1 Introduction 

In this paper we examine rigorously the phenomenon of quantum decoherence. This 
phenomenon is brought about by the interaction of a quantum system, called here "the 
system S" for short, with an environment, see e.g. [24, 36, 41] and the many references 
therein. Decoherence is reflected in the temporal decay of off-diagonal elements of the 
reduced density matrix of the system in a given basis. So far, this phenomenon has 
been analyzed rigorously only for explicitly solvable models, [15, 17, 24, 33, 36, 37, 35, 
40, 41]. In this paper we consider the decoherence phenomenon for quite general non- 
solvable models. Our analysis is based on the modern theory of resonances for quantum 
statistical systems as developed in [2, 6, 22, 9, 23, 21, 38, 31, 32] (see also the book 
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[20]), which is related to resonance theory in non-relativistic quantum electrodynamics 
([9, 5]). 

Let f) = f)s®l)R be the Hilbert space of the system interacting with the environment, 
also called the "reservoir", and let 

H = H s ® 1 R + l s ® H R + (1.1) 

be its Hamiltonian. Here, Hs and are the Hamiltonians of the system and the 
reservoir, respectively, and Xv is an interaction with a coupling constant A G R. We will 
omit trivial factors lis® and ®1r when there is no danger of confusion. In this paper 
we consider finite dimensional systems, the reservoirs are described by free massless 
quantum fields, and we take interactions of the form v = G (g> (f(g), where G is a 
self-adjoint matrix on f)s and (p{g) is the field operator, smoothed out with a coupling 
function (form factor). 

Consider possibly entangled initial states of the system and reservoir, where the 
reservoir is close to (a local perturbation of) an equilibrium state at some temperature 
T = 1/(3 > 0. (In the literature on decoherence, most often it is assumed that the 
initial states are product states, where the reservoir is in equilibrium, but our method 
works in the general case.) Let pt be the density matrix of the total system at time t. 
The reduced density matrix (of the system S) at time t is then formally given by 

Pt = Tr R/ 9 t , 

where TVr is the partial trace with respect to the reservoir degrees of freedom. For the 
sake of explicitness we describe here the case where the state of the reservoir is given 
by a well-defined density matrix on the Hilbert space In the next section we define 
the relevant notions for a more realistic reservoir, obtained for instance by taking a 
thermodynamic limit, or a continuous-mode limit. 

Let p((3, A) be the equilibrium state of the interacting system at temperature T = 
1/(3 and set p((3, A) := Tr^p((3, A). There are three possible scenarios for the asymptotic 
behaviour of the reduced density matrix, as t — > oo: 

0) Pt — ► Poo =p(P,\), 

(ii) Pt — ► Poo ^p(/3,A), 

(iii) ~p t does not converge. 

The first situation is generic while the last two are not, although they are of interest, 
e.g. for energy conserving, or quantum non-demolition interactions, characterized by 
[Hs,v] = 0, see [24] and Section 7.2. 

Decoherence is a basis-dependent notion. It is usually defined as the vanishing of the 
off-diagonal elements {p t ] m ,n, m ^ n m the limit t — > oo, in a chosen basis. Most often 
decoherence is defined w.r.t. the basis of eigenvectors of the system Hamiltonian Hs 
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(the energy basis, also called the computational basis for a quantum register), though 
other bases, such as the position basis for a particle in a scattering medium [24], are 
also used. 

Since p((3, A) is generically non-diagonal in the energy basis, the off-diagonal ele- 
ments of p t will not vanish in the generic case, as t — > oo. Thus, strictly speaking, de- 
coherence in this case should be defined as the decay (convergence) of the off-diagonals 
of p t to the corresponding off-diagonals of p(/3, A). The latter are of the order O(A) and 
in concrete applications often of the order 0(A 2 ). If these terms are neglected, then 
decoherence manifests itself as a process in which initially coherent superpositions of 
basis elements ipj become incoherent statistical mixtures, 

X^I^X^I — > ^PMj){^jV as t ^ oo. 
j,k j 

In particular, phase relations encoded in the Cj k, j 7^ k, disappear for large times. Of 
course, as t — > oo, the off-diagonal elements of p t vanish in a basis of eigenvectors of 
the asymptotic Hamiltonian #s,A,/3 := — \ hi/5(/3, A) = Hs + 0(A). We conjecture that 
this Hamiltonian absorbs the leading order correction to the non-interacting dynamics 
e~ ltHs due to the interaction with the reservoir. We discuss the role of i?s,A,/3 i n more 
detail elsewhere. We set H equal to one in what follows. 

In this paper we consider iV-dimensional quantum systems interacting in a quite 
general way with reservoirs of massless free quantum fields (photons, phonons or other 
massless excitations) . Let A be an arbitrary observable of the system (an operator on 
the system Hilbert space f)s) and set 

(A) t := Trs(p t A) = Tr s+R (p t (A 1 R )). (1.2) 

We show, under certain conditions on the interaction, that the ergodic averages 

(W)oo--=to* \ [ T (A) t dt (1.3) 

T->oo 1 J Q 

exist, i.e., that (A) t converges in the ergodic sense as t — > oo. Furthermore, we show 
that for any t > and for any < u>' < 

(A) t - ((A)) „ = ^2 e it£ Re{A) + O (A 2 e -5[ max 4im e }W/2]^j ^ (L4) 

where the complex numbers e are the eigenvalues of a certain explicitly given operator 
K(u'), lying in the strip {z G C | < Imz < uo' /2}. They have the expansions 

e = e W =e-A 2 5^ + 0(A 4 ), (1.5) 

where e G spec(^Ts <X> Is — Us ^ Hs) = spec(Hs) — spec(Hs) and the 5^ are the 
eigenvalues of a matrix A e , called a level-shift operator, acting on the eigenspace of 
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H$ lg — Is <8> Hs corresponding to the eigenvalue e (which is a subspace of f)s <S> t)s)- 
The level shift operators play a central role in the ergodic theory of open quantum 
systems, see e.g. [29]. We describe them in Section 4. The terms R S (A) in (1.4) 
are linear functionals of A and are given in terms of the initial state, po, and certain 
operators depending on the Hamiltonian H. They have the expansion 

R e {A)= Xm,nA m , n + 0(\ 2 ), 

(m,n)S/ e 

where I e is the collection of all pairs of indices such that e = E m — E n , the being 
the eigenvalues of Hg. Here, A rritn is the (m, n)-matrix element of the observable A in 
the energy basis of Hs, and the x m ^ n are coefficients depending on the initial state of 
the system (and on e, but not on A nor on A). 

Discussion of (1-4) ■ In the absence of interaction (A = 0) we have e = e € R, 
see (1.5). Depending on the interaction each resonance energy e may migrate into the 
upper complex plane, or it may stay on the real axis, as A / 0. The averages {A) t 
approach their ergodic means ((A)) if and only if Ime > for all £ ^ 0. In this case 
the convergence takes place on the time scale [Ime] -1 . Otherwise (A) t oscillates. A 

(s) 

sufficient condition for decay is that Im<5e < (and A small, see (1.5)). 

Remark about the error term in (1-4)- The restrictive condition Ime < lu'/2 < 
ir/ (3 in (1.4) which implies (3 < cA 2 , for some constant c, can be eliminated by using 
renormalization group methods as in [9, 34], and our results can be upgraded to hold 
uniformly in T = 1//3 — ► 0. This point will be addressed elsewhere. 

There are two kinds of processes which drive the decay: energy-exchange processes 
and energy preserving ones. The former are induced by interactions having nonvanish- 
ing probabilities for processes of absorption and emission of field quanta with energies 
corresponding to the Bohr frequencies of S (this is the "Fermi Golden Rule Condi- 
tion", [9, 18, 29, 31, 32]). Energy preserving interactions suppress such processes, 
allowing only for a phase change of the system during the evolution ( "phase damping" , 
[35, 12, 15, 17, 24, 33, 37]). 

Relation (1.4) gives a detailed picture of the dynamics of averages of observables. 
The resonance energies e and the functionals R £ can be calculated for concrete models, 
to arbitrary precision (in the sense of rigorous perturbation theory in A). See (1.9)- 
(1.11) for explicit expressions for the qubit, and the illustration below for an initially 
coherent superposition given by (1.15). In this paper we use relation (1.4) to discuss 
the processes of thermalization and decoherence. It would be interesting to apply the 
techniques developed here to the analysis of the transition from quantum behaviour to 
classical behaviour (see [11, 15]). 

We apply our results to a qubit, as well as to energy-preserving, or non-demolition 
interactions. They apply equally well to a register of arbitrarily many qubits. The case 
of energy-preserving interactions can be solved explicitly and serves as an illustrative 
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example as well as a starting point for a perturbation theory for interactions which are 
not energy-preserving, but for which the commutator [Hs,v] is small. 

Our results for the qubit can be summarized as follows. Consider a qubit coupled 
linearly to the field by the interaction 



v = 



a c 
c b 



<p(g), (i-6) 



where (f(g) is the Bose field operator, smeared out with a coupling function (form 
factor) g(k), k G R 3 , and the 2x2 coupling matrix (representing the coupling operator 
in the energy eigenbasis) is hermitian. The operator (1.6) - or a sum of such terms, for 
which our technique works equally well - is the most general coupling which is linear 
in field operators. We refer to Remark 14 below for the link between (1.6) and the 
spin-boson model. 

Note that the form-factor g contains an ultra-violet cut-off which introduces a time- 
scale tjjv- This time scale depends on the physical system in question. We can think 
of it as coming from some frequency-cutoff determined by a characteristic length scale 
beyond which the interaction decreases rapidly. For instance, for a phonon field tjjv 
is naturally identified with the inverse of the Debye frequency. We assume tjjv to be 
much smaller than the time scales considered here. 

A key role in the decoherence analysis is played by the infrared behaviour of form 
factors g G L 2 (IR 3 , d 3 k). We characterize this behavior by the unique p > —1/2 satis- 
fying 

< lim = C < oo. (1.7) 

IfcHo \k\P V ; 

The power p depends on the physical model considered, e.g. for quantum-optical 
systems, p = 1/2, and for the quantized electromagnetic field, p = —1/2. 

Decoherence of models with interaction (1.6) with c = is considered in [12, 15, 17, 
24, 33, 35, 37, 39] (see also Section 7.2). This is the situation of a non-demolition (energy 
conserving) interaction, where v commutes with the Hamiltonian H$ and consequently 
energy-exchange processes are suppressed. The resulting decoherence is called phase- 
decoherence. A particular model of phase-decoherence is obtained by the so-called 
position-position coupling, where the matrix in the interaction (1.6) is the Pauli matrix 
a z [12, 17, 35, 39]. On the other hand, energy-exchange processes, responsible for 
driving the system to equilibrium, have a probability proportional to \c\ 2n , for some 
n > 1 (and a, b do not enter) [1, 9, 18, 22, 29, 31, 32]. Thus the property c / is 
important for thermalization (return to equilibrium). 

We express the energy-exchange effectiveness in terms of the function 

«"» = 5?; l d ' Wh (^) lsik)f \m- e ,r +c - (L8> 

where r] > represents the energy at which processes between the qubit and the 
reservoir take place. Let A = E2 — E\ > be the energy gap of the qubit. In 
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works on convergence to equilibrium it is usually assumed that |c| 2 £(A) > 0. This 
condition is called the "Fermi Golden Rule Condition". It means that the interaction 
induces second-order (A 2 ) energy exchanging processes at the Bohr frequency of the 
qubit (emission and absorption of reservoir quanta). The condition c / is actually 
necessary for thermalization while £(A) > is not (higher order processes can drive 
the system to equilibrium). Observe that £(A) converges to a fixed function, Co (A), as 
T — > 0, and £(A) increases exponentially as T — > oo. The expression for the decoherence 
involves also £(0) (see below). 

In this paper we describe the dynamics, and in particular the decoherence prop- 
erties, of systems which exhibit both thermalization and (phase) decoherence. See the 
the discussion after (1.7) for a comparison of the two effects, and how they relate to 
the coefficients a, b, c and the coupling function g in (1.6). 

Let the initial density matrix, pt=o, be of the form ~p <g> pu p. (Our method does 
not require the initial state to be a product, see Remark 5 below.) Denote by Pm,n 

the 

rank-one operator represented in the energy basis by the 2x2 matrix whose entries 
are zero, except the (n, m) entry which is one. We show that for t > 

fo]i,i - «Pi,i»oo = e ife0(A) [Co + 0(A 2 )] (1.9) 

+ e iteA(A) 0(A 2 ) + e ite - A(A) 0(A 2 ) + 0(X 2 e~ tuJ ' /2 ) 

and 

[ft]i, 2 -((pi,2)) 00 = e ite -W[C A + 0(A 2 )] (1.10) 

+ e ite0(A) O(A 2 ) + e ite - A(A) 0(A 2 ) + 0(A 2 e"^ 72 ). 

Here, Co,Ca are explicit constants depending on the initial condition p , but not on 
A, and the resonance energies e have the expansions 

e (A) = iAV|c| 2 £(A) + 0(A 4 ) 

e A (A) = A + A 2 i? + |AV [| c | 2 £(A) + (6-a) 2 £(0)] + 0(A 4 ) (1.11) 
£-a(A) = -£a(A) 

with the real number 

R=i(b 2 -a 2 )(g,u>- 1 g) 
+I|c| 2 P.V. f u 2 \g(\u\,a)\ 2 cotJ^)' 
The error terms in (1.9), (1.10) and (1.11) satisfy, for small A, 



0(X 2 ) , , 

< C and sup 



A 2 



0(A 2 e-*-'/2) 



\2 e -tw'/2 



< C. 



Remarks. 1) To our knowledge this is the first time that formulas (1.9)-(1.11) 
are presented for models which are not explicitly solvable. Results for exactly solvable 
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models (non-demolition interactions) are given, among others, in [35, 37, 33]. See [24] 
for an overview of the subject. 

2) Relations (1.9)-(1.11) are valid for all values oft > 0, and the remainder terms 
0(A 2 ) are uniform in t > 0. In particular, we do not require that A — > as t — > oo (van 
Hove limit). 

3) Even if the initial density matrix, pt=o, is a product of the system and reservoir 
density matrices, the density matrix, pt, at any subsequent moment of time t > is 
not of the product form. In other words, the evolution creates the system-reservoir 
entanglement. 

4) The corresponding expressions for the matrix elements [p t ]2,2 and [p t ]2,i are 
obtained from the relations [pth,2 = 1— [Pt]i,i (conservation of unit trace) and [p t ]2,i = 
\Pt]\,2 (hermiticity of p t ). 

5) If the qubit is initially in one of the logic pure states ~p = {p^iPjl, where 
H s <pj = Ejipj, j = 1, 2, then we have C A = 0, and C = e^ A / 2 (e^ A + 1)~ 3 / 2 for j = 1 
and C = e^ A (e^ A + l)" 3 / 2 for j = 2, see at the end of Section 7.1. 

6) We develop a formula for (A) t — ((A)) ^ for all observables A of any iV-level 
system S in Section 4. 

7) If the system has the property of return to equilibrium, i.e., if £(A) > 0, then 

((Pn,m)) = [Poo]m,n = 0~m,n ^ ^ e ~f3H s j + )' 

We thus recover the Gibbs law in the long time limit followed by the weak coupling 
limit. A similar observation is found in the context of the quantum Langevin equation 
in [10]. 

8) If po is an arbitrary initial density matrix on Tt$ <g> TL^ (i.e., not necessarily of 
product form), then the method developed in Section 4 yields the following result: For 
any rj > there are constants Co, Ca, depending on r] and po but not on A, s.t. 

foh,i-«Pi,i»oo = e ifeo(A) [Co + O(X)} (1.12) 

+ e ife AWo(A) + e ite - A(A) 0(A) + 0(A e"^' /2 ) + 0(rj) 

and 

[p t ]i,2 ~ ((pi,2))oo = e ite -W [C A + 0(A)] (1.13) 

+ e iteo(A) 0(A) + e it£ -^0(X) + 0(Ae^^' /2 ) + 0(rj), 

where 0(r]) is uniform in t, and where the resonance energies are given by (1.11). 
Furthermore, all remainder terms depend on rj, in general. 

9) Equations (1.9), (1.10) and (1.11) define the decoherence time scale, td = 
[ImeA(A)] -1 , and the thermalization time scale, tt = [Imeo(A)] -1 . We should com- 
pare td with the decoherence time scales in real systems and with computational time 
scales. The former depends on the physical realization of the qubit and its environment. 
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It can vary from 10 4 s for nuclear spins in paramagnetic atoms to 10~ 12 s for electron- 
hole excitations in bulk semiconductors (see e.g. [16]). ReeA(A) — A = X 2 R + 0(A 4 ) 
gives the radiative energy shifts. 

10) To second order in A, the imaginary part of ea is increased by a term oc 
(b — a) 2 £(0) only if p = —1/2, where p is defined in (1.7). For p > —1/2 we have 
£(0) = and that contribution vanishes. For p < — 1/2 we have £(0) = oo. 

11) £(A) and R contain purely quantum, vacuum fluctuation terms, as well as 
thermal ones, while £(0) is determined entirely by thermal fluctuations. £(A) and £(0) 
are increasing in T, and, as T — > 0, £(0) is linear in T (p = — 1/2) and £(A) converges 
to a fixed nonzero value. The decoherence rate thus increases for decreasing T, and 
it approaches a finite value as T — > 0, for c / 0. Our proofs work for arbitrarily 
small, fixed temperatures. There is strong evidence that the results above remain valid 
for T — > (see the remark about the error term in (1.4) above, and also [9, 34]). A 
discussion of the decoherence function in terms of the temperature for the explicitly 
solvable case, c = 0, is given in [35]. 

12) To second order in pertrubation, the ratio of the thermalization versus deco- 
herence rate is tt/td = |[1 + (^p) 2 |^]- For tt/td < 1, the populations converge to 
their limiting values faster than the off-digaonal matrix elements, as t — > oo (coherence 
persists beyond thermalization of the population). For tt/td > 1, the off-diagonal el- 
ements converge faster. If the interaction matrix is diagonal (c = 0), then tt/td = oo, 
if it is off-diagonal (or if a = b), then tt/td = 1/2. 

13) For energy-conserving interactions, c = 0, it follows that full decoherence occurs 
if and only if b ^ a and £(0) > 0. If either of these conditions are not satisfied then the 
off-diagonal matrix elements are purely oscillatory (while the populations are constant). 
We analyze energy-conserving interactions in Section 7.2. 

14) In the ubiquitous spin-boson model [28], obtained as a two-state truncation 
of a double-well system or an atom, interacting with a Bose field, the Hamiltonian 
is given by (1.1) with H$ = — ^Aoa x + \to z and v = a z f{g)- Here, a x , a z are 
Pauli spin matrices, e is the "bias" of the asymmetric double well, and Ao is the "bare 
tunneling matrix element". In the canonical basis, whose vectors represent the states 
of the system localized in the left and the right well, Hs has the representation 



b 2 



e -A 
-A -e 



(1.14) 



The diagonalization of H s yields H s = diag(E + ,E^), where E± = ±\\/ e 2 + A^. The 
operator v = <j z ip(g) is represented in the basis diagonalizing Hs as (1.6), with 
a = -b = -(^ + l)- x /2 an d c = *(^ + l)-V2. 

Illustration. Let us discuss the decoherence and thermalization properties in the 
case where S is initially given by a coherent superposition in the energy basis. For sake 
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of explicitness we take 

Po = h \\ • (i-is) 

We obtain the following expressions for the dynamics of matrix elements, for all t > 0: 

e -(3E m ( — l) m / 3A\ 

[ptUm = + tanh e lte °( A ) + i? m , m (A, t), m = 1,2, 

[Pt]i,2 = ie ite --W + J R li2 (A,t), 

where the numbers e are given in (1.11). The remainder terms satisfy |i? mi „(A,t)| < 
CA 2 , uniformly in t > 0, and they can be decomposed into a sum of a constant and a 
decaying part, 

e -/3E m f 

where |i?^„(A,t)| = 0(X 2 e~' ft ), with 7 = min{Imeo, Ime-|-A}. These relations show 
in particular that, to second order in A, convergence of the populations to the equi- 
librium values (Gibbs law), and decoherence occur exponentially fast, with rates 
tt = [Imeo(A)] -1 and td = [Im£A(A)] _1 , respectively. (If either of these imaginary 
parts vanishes then the corresponding process does not take place, of course.) In 
particular, coherence of the initial state stays preserved on time scales of the order 
A -2 [|c| 2 £(A) + (b — a) 2 ^(0)] _1 , c.f. (1.11). We show how to arrive at the above expres- 
sions at the end of Section 7.1. 

The method we use in this work yields an error estimate in (1.4) which is not 
uniform in T = 1/(3 — > 0. However, our result can be upgraded to a uniform estimate 
by employing spectral renormalization group methods as developed in [9, 31, 32, 34]. 
This will be addressed elsewhere. 

As mentioned above, we prove equation (1.4) using quantum statistical resonance 
theory, which is based on spectral deformation techniques. To keep the exposition as 
simple as possible we use in this paper the simplest and most restrictive version of this 
method, namely, the translation deformation. This produces weaker results (like non- 
uniformity in temperature of the error estimate in (1.4), as mentioned above) than can 
be obtained by more refined techniques. Dilation deformation [9] and a combination 
of dilation and translation deformation [31, 32] in conjunction with spectral renormal- 
ization group methods weakens the restrictions on the class of treatable interactions 
considerably and strengthens the results, but at the price of a much more involved 
mathematical machinery. 

Although there is a subtle mathematical theory behind our techniques, we believe 
that on a formal level, they are simpler and more powerful than the standard path- 
integral methods ([14, 3, 25, 26]). The rigorous treatment based on master equations 
and Lindblad generators [27] seems to be more difficult. While the path-integral and 
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master equation approaches are intrinsically time-dependent, the resonance theory is 
formulated as a stationary eigenvalue problem for some (albeit non-self-adjoint) oper- 
ators ^ e ,/3,x acting on subspaces of two copies of the system Hilbert space, f)s <8> f)s- m 
particular, the complex numbers e are eigenvalues of these operators. See the remark 
at the end of Section 4. 

This paper is organized as follows. In Section 2 we introduce the model and in 
Section 3 we present our main result for general N— level systems and for the qubit. 
We develop a general dynamical resonance theory for systems at positive temperatures 
and densities in Sections 4-6 (our general theorem stating (1.4) is proven in Section 4 
and the relation (1.5) is shown in Section 5). We give two applications of this theory 
in Section 7: the first one yields a proof of the results for the qubit mentioned above, 
the second one illustrates our method on an explicitly solvable non-demolition model. 
Appendices A and B contain explicit formulas for some quantities which are important 
for our resonance theory. In Appendices C and D we outline the perturbation theory of 
equilibrium states, and we give the proofs of several propositions of previous sections, 
including the short calculation of the explicit solution of a non-demolition model. 

2 Open quantum systems 

The system S. Let f)s be the space of pure states of a quantum system S whose 
dynamics is generated by a Hamiltonian H$. In applications we will take f)s = C N and 
Hs = diag(-Ei, . . . , En), but we give a discussion of more general systems S with finitely 
or infinitely many (discrete) levels. This discussion is straightforward and instructive, 
we believe. The evolution of a density matrix ps on fjg is 

t ' ^ e- itHs p s j tHs , t G M. (2.1) 

The Gibbs state at temperature T = 1/(3 > is ps,/3 = ^ 1 e _/3Hs , with the nor- 
malization constant Zp = Trse~^ s . (If dimfjs < oo then Zp < oo for any H§. In 
the infinite-dimensional case we assume that is trace-class.) Observables of S are 
operators on rjs, they form the algebra of all bounded operators H(f)s)- 

One can represent any state on S, mixed or pure, by a single vector in the Hilbert 
space Hs = f)s ® f)s- This is the so-called Gelfand-Naimark-Segal representation of 
states. To see how this works take an arbitrary density matrix ps on fjs, and write it 
in its diagonalizing basis as 

PS = Z Zpj\^)(^\- (2-2) 

3 

The corresponding state is defined by A ^ Trg(ps^4) for any observable A G £>(hs). Let 
us denote by L 2 (f)s) the set of all Hilbert-Schmidt operators, i.e., A G £ 2 (f)s) if an d 
only if Tts(A*A) is finite. (In the case dimf)s < oo we have -L 2 (rjs) = #(f)s)0 £ 2 (f)s) is 
a Hilbert space with the scalar product (A,B) 2 = Trs(^4*.B). Since the density matrix 



10 



ps is a positive trace class operator, its square is a Hilbert-Schmidt operator, and due 
to the cyclicity of the trace we have 

Tr s (psA) = {y/te,Ay/pd 2 . (2.3) 

This shows that the density matrix ps on f)g is represented by the vector y/ps G L 2 (f)s)- 
Instead of working with the state space -L 2 (f)s) it is often convenient to switch to 
Hs = f)S <8> f)S- This is done via the correspondence 

T ■ \<p)(ip\ — ► f ® Cip, (2.4) 

which extends by linearity to an isometric isomorphism between -L 2 (f)s) and Hs- Here, 
C is a map on f)g that is chosen to be antilinear (conjugate linear) since T should be 
linear but ip i— > (^| is antilinear. Furthermore, to make T an isometry (norm preserving) 
C has to be antiunitary, meaning that it satisfies 

(C X i,CX2) =Tx^X2) (2.5) 

for all xi,2 G f)s- Relation (2.5) implies that C is bijective: injectivity follows from 
ll^xll = llxll and surjectivity follows from the fact that C 2 is unitary, so RanC 5 
RanC 2 = fj s . 

In what follows we choose C to be the operator that takes the complex conjugate 
of coordinates w.r.t. the basis of f)s in which H$ is diagonal. In this case C is an 
involution, C 2 = C. 

Since T is an isomerty we have (A, A) 2 = ||TA|| 2 = (TA,TA), where the norm and 
the scalar product on the r.h.s. are those of Hs- This implies that (A, B) 2 = (TA, TB), 
by the polarization identity. In particular, (2.3) gives 

Tr s (ps^4) = (T^/ps, TA^) . (2.6) 
The vector in Hs representing the state (2.2) is thus 

y = T^ps- = Y,Vp]i>j® c i>j, ( 2 - 7 ) 

3 

and the Gibbs state corresponds to 

Qs,(3 = Z~ 1/2 e- fiE ' /2 <Pj ® <Pj, (2.8) 

3 

where (fj is the eigenvector of Ho, corresponding to the energy Ej. 
Due to (2.6) and since TA^/ps = (A <g) 1)^ we have 

Trs(ps^) = <*,(A®1)*). (2.9) 
The algebra of observables of S is given by £>(f)s) <8> 1 when viewed as operators on Hs- 
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Remark. It appears that in the representation of the Gibbs state the temperature 
dependence (parameter (3) is entirely concentrated on the vector fisfl) (2-8), and the 
represented observables A®! are independent of (5. We may transfer the /3-dependence 
entirely (or partly) to the observables by a change of basis. For example, let be any 
fixed vector on Hs and take any unitary Up with the property UpQs,p = Then we 
have Trs(ps-4) = (*', Ap *'), where Ap = Up(A <g> t)U* p . 

As we see from (2.1), (2.2) and (2.7), the evolution of * is given by 

t^J2VPj e~ HH *ipj ®Ce- itHs ipj = e~ itLs ^, (2.10) 

3 

where 

L s = H s ®t-t® H s (2.11) 
is called the (standard) Liouville operator. It satisfies the relation 

L s n St p = 0. (2.12) 

In the Heisenberg picture, a system observable A evolves according to 

t ' — ^ e itLs (A®t)e- itLs . (2.13) 

Remark on the choice of the Liouville operator. The map T, (2.4), is only a 
particular choice of an isometric isomorphism. We may more generally put 

T : \cp)(ijj\ — ► U<p® VCip, 

where U and V are arbitrary unitary operators on f)s- We may even define a time- 
dependent isometric isomorphism Tj , by taking time-dependent families of unitaries Ut 
and Vf. For instance, the special choice [7 = 1 and V = Vt yields 

T t Ae- itHs ^ps~j tHs = (A (g) t)(e~ itHs ® V^e^ 3 )^, 

with ^ given in (2.7). Therefore, we may equally well define the evolution of ^ by 
t i > = (e _i * Hs (8) Fte 1 *^)^. As an example, take V t = e~ itHs , then *(i) = 

(e-' ltHs (g> 1)*. It is easy to see that the only choice of Vt for which the Gibbs vector 
fig p, (2.8), is invariant under the evolution, i.e. for which S7g pit) = ^s,/3 f° r an ^ is 
given by Vt = t. This convenient choice results in (2.10), (2.11) and (2.12). 

The reservoir R. We introduce now a second quantum system, a reservoir R 
in a state of thermal equilibrium at a temperature T = l/{3 > 0. A typical reservoir 
is a very large quantum system, say a quantum gas with a given particle density (or 
a radiation field). In order to suppress recurrences one may consider the idealized 
situation of an infinitely extended reservoir giving rise to truly irreversible processes 
(like thermalization or decoherence) . In physical experiments reservoirs have of course 
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finite size, but the idealized limit is a good approximation for physically realistic times 
that are not exceedingly large, see e.g. [19] for a discussion of this point. 

The reservoir R we consider is the infinitely extended gas of massless relativistic 
Bosons (of photons or phonons, for example) at positive temperature and positive den- 
sity. Our approach applies also to fermionic reservoirs, and in fact becomes technically 
much simpler in that case. The state of R is obtained by performing a thermodynamic 
limit of finite-volume equilibrium (Gibbs) states with fixed temperature and density. 
We refer to [30] for a detailed exposition of this matter in textbook-style. Just as in 
the case of the system S one can represent the equilibrium state by a single vector in 
a suitable Hilbert space. This is the so-called Araki- Woods representation, [4, 30], in 
which the Hilbert space is given by 

n R = F(L 2 {M. 3 , d 3 fc)) ® F(L 2 (R 3 , d 3 k)), (2.14) 

where J r (L 2 (IR 3 , d 3 /c)) is the bosonic Fock space over the one-particle space of wave 
functions L 2 (R 3 ,d 3 /c) in momentum representation. The usual bosonic creation and 
annihilation operators a*(k) and a(k), k £ IR 3 , are represented in TC R as the thermal 
creation- and annihilation operators, 

a(k) i-> ap(k) = + fip a(k) <g> t + JJTp 1 <g> a*(k), (2.15) 
a*(k)^a*p{k) = y/1 + [ip a*(k) <g> 11 + y/JTp 11 <g> a(k), (2.16) 

where fip is Planck's momentum density distribution for black-body radiation, 

M*) = ^TT ( 2 - 17 ) 

(we take the Bose gas in a phase without a Bose-Einstein condensate). The represen- 
tation (2.15), (2.16) is the equivalent to the representation of an observable A G S(f)s) 
by A <S> l £ B(Hs) in the case of the system S, c.f. the remark after (2.9). The map 
defined by (2.15), (2.16) is called the Araki- Woods representation map. 

Denote by Q the vacuum vector of ^ r (L 2 (M 3 , d 3 fc)). It is easily seen that the vector 

n R) p = n<g>n (2.18) 

represents the equilibrium state and 

(^R,/3, a*p(k)ap(l)n Rt p) = 5(k - l)/ip(k), 

where 5 is the Dirac delta function. 

The dynamics of a density matrix p R of the reservoir (acting on Hr) is given by 

t ^ e- itLR p R e itLR , (2.19) 

where 

£r = H R (8> 11 - H <g> H R , (2.20) 
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with 



H R = f \k\a*(k)a(k)d 3 k. (2.21) 



The self-adjoint operator Lr is called the standard Liouville operator for R, it 
satisfies (compare with (2.12)) 

£r^r 5/ 3 = 0. (2.22) 
The smoothed-out thermal creation and annihilation operators are defined as 

a*Jh) = [ h(k)a* g (k)d 3 k, (2.23) 

JWt? 

ap(h) = [ h(k)a p (k)d 3 k, (2.24) 

where h G L 2 (M. 3 ,dk 3 ) is a wave function of a single Boson and where the ag(/c) and 
a*p(k) are given in (2.15), (2.16). 

It is not hard to check explicitly that Lr, (2.20), implements the Heisenberg dy- 
namics of observables, given by the Bogoliubov transformation 

1 1 ^ e itL *a*{h) e~ itLR = aj( e^/i). 

The total system S + R. The joint system S + R is described by the Hilbert 

space 

H = H S ®H R , (2.25) 
and the non-interacting dynamics of a density matrix p on 7is <8> Wr is 

1 1 — ^ e - itL Ve itLo , (2.26) 

with 

L = L s + L R . (2.27) 

The state ps,/3 ® Pr,/3 = |^s,/3)(^S,/?| ® |^R,/3)(^R,/3| is an equilibrium state w.r.t. the 
non- interacting dynamics. 

The coupling between S and R is specified by an interaction operator V which is 
an observable of the joint system and is related to the interaction operator v given in 
(1.6). The full dynamics of a density matrix p on 7is (8> Wr, 

e - itL *pe itL \ (2.28) 

is generated by 

L x = L + XV, (2.29) 

where A is a coupling constant. An important class of interactions is given by coupling 
operators of the form 

v = G®ip(g), (2.30) 



14 



see also (1.1), where G is a self-adjoint operator on f)g and where 

i p(g) = -L(a*(g)+a(g)) (2.31) 

is the bosonic field operator on JF(L 2 (M 3 , d 3 /c)), smeared out with a form factor g G 
L 2 (R 3 , d 3 k). On the positive temperature Hilbert space Hs <8> Hr this interaction 
Hamiltonian leads to the operator 

V = G®l s ®<pp(g), (2-32) 

where 

tp /3 (g) = -j=(a* p (g) + a p (g)) (2.33) 

is the smoothed out thermal field operator, (2.23), (2.24). More generally, the interac- 
tion is specified by an operator V £ Wl = Wis ® SDTr, or a V which is unbounded, but 
affiliated with Wl, 1 like (2.32). 

For a wide class of interactions V, the interacting dynamics admits an equilibrium 
state p(/3, A) (at least for small A), but we do not limit our discussion to such operators 
V. 

The reduced density matrix ~p t of the open system S is obtained by tracing out the 
degrees of freedom of the reservoir, it is defined by 

Tr s (p t A) = Trg+R ( e^Vo e ltL "(A 1 R )) , 

for all A e B(t)s), where po is the initial state of S. It is customary to take po = 
Po ® PR,/9) where ~p is the initial state of S and pr^ is the reservoir equilibrium state, 
but our analysis works equally well for initially entangled (non-product) states which 
are arbitrary local perturbations of Pr q- 

Let {fj}j>i be an orthonormal basis of f)s diagonalizing Hs- The reduced density 
matrix has matrix elements 

[Pt\m,n ■= {<Pm,Pt<Pn} = T^siPtP^m) 

where 

Pn,m = \tPn)(<Pm\- (2-34) 

We say that the system S + R has the property of return to equilibrium iff 
lim Trg+R^oe^Me-^) = Tr s+R (p(/3, A)M), 

t^oo 

for all observables (of the joint system) M and for all initial density matrices po on H. 
The large time limit of the reduced density matrix of such a system is given by 

p^ := lim p t = Tv R {p{/3, A)) = p s ,p + 0(A), 

t — >oo 

1 meaning that it commutes with all bounded operators belonging the the commutant algebra QJl' 
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since, by perturbation theory of equilibrium states, p({3, A) = ps,/3 ® PR,p + 0(A). The 
leading term of (for small coupling) is just the Gibbs state of S. In this sense, S 
undergoes the process of thermalization. 

The system S is said to exhibit (full) decoherence (in the energy basis) if the off- 
diagonal matrix elements of the reduced density matrix vanish in the limit of large 
times, 

lim [p t ] m ,n = 

whenever m / n. 

It has to be pointed out that is not diagonal in the energy basis in general, but 
only its leading part is. Therefore, thermalization prevents full decoherence since the 
final state of the thermalized system S has still non-vanishing off-diagonal elements, 
which are of the order 0(A). In fact, there is creation of non-vanishing off- diagonals by 
thermalization since even if S is initially an incoherent superposition of basis elements 
of a given basis {Xj}j>ii 

Po = ^2a m \xm)(Xm\, 

m 

it will become a coherent superposition since in general, 

(Xm,PooXn) + 0. 

We illustrate this for a qubit in Theorem 3.3. 

Full decoherence can occur if thermalization processes are excluded. For instance, if 
one suppresses energy exchanges between S and R (i.e., if V commutes with Hs) then, 
even if an equilibrium state exists, the system does not converge to it as t — > oo. So 
there is no thermalization but full decoherence, i.e. [Poo]m,n = 0, m ^ n, can be shown 
to occur in various models. (That there is no thermalization is readily seen since the 
populations are time- independent). 

However, it is well known [35, 33] that even if energy exchange processes are sup- 
pressed, full decoherence may not take place if the infrared modes of the reservoir are 
only weakly coupled to the system, see Section 7.2. 

3 Decoherence and thermalization of a general iV-level 
system and of a qubit 

We analyze an iV-level system coupled to a thermal environment modeled by an in- 
finitely extended free massless relativistic Bose field at temperature T > 0. Under 
suitable conditions on the interaction the system has the property of return to equi- 
librium. Our goal is to compare the thermalization and decoherence processes and, in 
particular, the speed of convergence of the diagonal and off-diagonal matrix elements 
of the reduced density matrix. 



16 



The state space of an iV-level system is given by hs = C , and its Hamiltonian is 

H s = diag(E 1 ,...,E N ). (3.1) 

We couple the iV-level system to the reservoir through an operator v = G®tp(g), where 
G is a hermitian N x N matrix, and the field operator is as in (2.31), with a coupling 
function g{k), k £ M 3 . Let g(k) = g(r, a), where (r, a) 6 M+ x S 2 . We make a regularity 
assumption on g. Fix any phase (f> € M and define 

/■ \ / ^ i ,1/2 f 9{ u ,°~) if tt > 0, ,„ „s 

9 ^ a) ■= ii^-u m 1/2 { if « < o! (3 - 2) 

where n£l and a £ S* 2 . The phase is a parameter which can be chosen appropriately 
to satisfy the following condition for a given coupling function g, see [18]. 

(A) We assume that the map uj ^ gp(u + uj,a) has an analytic extension to a 
complex neighbourhood {\z\ < uj'} of the origin, as a map from C to L 2 (M 3 , d 3 k). 

This condition ensures that the simplest version of a dynamical resonance theory 
- the one using complex deformations - can be implemented in a straightforward way. 
Examples of g satisfying (A) are given by g(r, a) = r p e~ r ' m g\(a), where p = — 1/2 + n, 
n = 0, 1, . . ., m = 1, 2, and gi(a) = e 1< >'g 1 (a). 

The technical simplicity of the complex translation method comes at a price. On 
the one hand, it limits the class of admissible functions g(k), which have to behave 
appropriately in the infra-red regime so that the parts of (3.2) fit nicely together at 
u = 0, to allow for an analytic continuation. On the other hand, the square root in 
(3.2) must be analytic as well, which implies the condition uj 1 < 2tt/ (5. 

We now state our main result on the dynamics. Let tpQ <S> ^r,/3 be the vector in 
TCs <8> 7~Cr representing the density matrix at time t = 0. Let B be the unique operator 
in il s ® B(C N ) satisfying 

Btt s ,p = ^o- 

The existence of such a B is not hard to verify, see also the end of Section 7.1 for 
concrete examples. Define the vector 

Theorem 3.1 (Dynamical resonance theory) Assume condition (A) with a fixed 
uj' satisfying < uj' < 2n/f3. There is a constant Co s.t. if |A| < cq/P then the limit 
((^4)) 0O , (1-3), exists for all observables A G £>(f)s)- Moreover, for all such A and for 
all t > we have 

(A) t -((A)) 00 = £ j^ s) ((B*^ )®n R> p,Qi*\A®l s )np, ) 

"(c) 

+ EE ^ ((B*4>o)®n Rt p,Ql s \A® l s )%o) 

e^0 s=l 

+0(A 2 e _ 5[ maXe . s i Im } £ c s) +w72]^ ^33-) 
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Here, the are given by (1-5), 1 < v(e) < mult(e) counts the splitting of the eigen- 
value e into distinct resonance energies e*i\ and the are (non-orthogonal) finite- 
rank projections. 

Remarks. 1. One can explicitly expand in powers of A both the resonance energies 
e*i^ and the resonance projections Q^e \ see Sections 5 and 6. 

2. Our techniques are not restricted to unentangled initial states (see also (1.12), 
(1.13)). 

Return to equilibrium for TV-level systems with a coupling (2.32) to the Bose field 
has been shown by several authors over the last few years, provided the coupling is 
sufficiently effective, and under regularity conditions much weaker than (A), see for 
instance [22, 9, 18] and references therein. For the sake of notational simplicity we 
will restrict our attention to the qubit {N = 2) in the presentation of the remaining 
results in this section. The Hamiltonian of the qubit is H$ = diag(.Ei, E2), we set 
A := E2 — Ei > 0. The operator G in (2.32) is represented in the energy-basis by the 
matrix 

r a c 
c b 



G = 
where a,6el. 



(3.4) 



Theorem 3.2 (Thermalization, return to equilibrium for the qubit) Assume 
condition (A), and that |c| 2 £(A) > 0. There is a constant cq s.t. if < |A| < co/(3 
then the qubit coupled to the reservoir of thermal Bosons has the property of return to 
equilibrium. 

The condition |c| 2 £(A) > is often called the "Fermi Golden Rule Condition". It 
means that the interaction induces second-order (A 2 ) energy exchanging processes at 
the Bohr frequency of the qubit (emission and absorption of reservoir quanta). The 
condition c / is actually necessary for thermalization while £ (A) > is not (higher 
order processes can drive the system to equilibrium). 

Theorem 3.3 (Creation of off-diagonals by thermalization) 

Under the assumption of Theorem 3.2 the off-diagonal elements of the reduced den- 
sity matrix are given in the limit t — > 00, irrespectively of the initial density matrix po, 
by 

[Pji,2 = ^- (9, (ae~^B 1 + be-^B 2 )g) + 0(A 4 ), 
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where Zs t p 



Trg(e P H s) is the partition function ofS, and where 



Bi 




and 



B 2 



( e -/9w/2 _ e /?A/2^2 + e /3A/2 _ l 



lo(uj + A) . 



e ^_ e /3A/2 (e/ 3A/2_ 1) _ 1 + 



e/? A / 2 (e^ A / 2 -l) 



+ 



uj(lo - A) 
1 + H 



uA 



A(u + A) 



A(u - A) 



Here to = \k\, A = E2 — E\ > is the energy gap of Ho,, and p = p,p is given in (2.17). 

This theorem follows from an expansion of the equilibrium state in powers of A. We 
give an outline of the proof in Appendix C. 

Remarks. 1) Even if the system S starts out in an incoherent superposition of 
vectors from the energy basis, Theorem 3.3 shows that in the large-time limit, coherence 
of order A 2 is built up in that basis. 

2) If the interaction matrix G is diagonal, i.e. if c = 0, then the 0(A 2 ) terms in 
the off-diagonals of vanish. The same happens if a = b = 0, i.e., if G is purely 
off-diagonal. 

Recall the notation (A) t , ((A))^ and p n>m , (1.2), (1.3) and (2.34). 

Theorem 3.4 (Convergence of matrix elements) Assume (A) for a fixed to' > 
satisfying < u' < 2n/f3. There is a constant cq s.t. the following statements hold for 
|A| < co//3. For any initial density matrix of the formJ/Q® p^p and any observable A of 
S the limit ((^4)) 00 exists, and statements (1.9)-(1.11) about the asymptotic behaviour 
of [p t ]m,n are true. 

Let us again point out that our method works for any entagled initial density matrix 
p on Hs <8> Hr, c.f. (1.12), (1.13). 

4 Dynamical resonance theory for systems at positive 
temperatures and positive densities and proof of Theo- 
rem 3.1 

The goal of this section is to arrive at relation (3.3) expressing the asymptotic behaviour 
of the averages of a system observable A as t — > 00. 
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The state of S, given by the initial density matrix p , is represented by the vector 
state A i-> (tpo, (A <g) is) tpo), for a ipo £ = fys^S; see (2.7). For instance, -^o = ^s,/3 
if S is initially in equilibrium, c.f. (2.8). 

Let the initial state of the total system S + R be the product state ~p ® pR t p, which 
is represented on Hs <8> Wr by the vector 

® ^r,/3 e H s ® Wr, (4.1) 

where Or^ is given in (2.18). An observable A G 97T of the system S + R evolves 
according to the Heisenberg evolution 

where L,\ is given in (2.29), (2.32). The Schrddinger dynamics of a vector ip G <8>Hr 
determining the state ■ "0) on 971 is given by i i— > Vt = e~ ltLx ip. The average of an 
observable A of S at time i is thus 

(A) t = <^ ® fiR,/3, e itLA (4 ® Is ® 1r) e- itL ^ Vo ® ^R,/3> ■ (4.2) 

It is reasonable to just write A instead of A Is ® 1r in this section. Our goal is to 
examine (A) t as a function of t. 

It is not hard to see that ^s,/3 nas the following property of separability: given any 
vector ipQ £ H.s <g) Ws we can find a unique operator B G Is ® •B(f)s) with the property 
that 

= £%/3. (4.3) 

One can solve this equation for the matrix elements of B, e.g. by using (2.8). See also 
the end of Section 7.1 for concrete examples. 

An entangled initial state of the system, given by a density matrix po , is realized on 
the Hilbert space Hs ® Wr by a vector which is not of the product (4.1). However, 
one can show that any ^ can be approximated arbitrarily well by vectors of the form 
BQpfi, for some B in the commutant of 97t, and where 

ty9,0 = %/3 ® ^R,/3- (4.4) 

(The commutant of 97t is the von Neumann algebra consisting of all operators on 
TCs ® Wr which commute with all operators of 971.) The following arguments can then 
be carried out in a similar fashion. This is why our method works for arbitrary initial 
states, not only for Feynman- Vernon type, disentangled initial states. 

Since e ltLA 97te~ ltiA = 97? (invariance of the algebra of observables) , the operator B 
(belonging to the commutant of 971,) commutes with e ltLx (A® is <g> 1r) e~ ltLx . We can 
thus use (4.3) in (4.2) and commute B to the left, with the result 

(A) t = ((£*Vo) ® ^r,/3, e itLx A e- itLA %o> • (4.5) 
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The following formula, derived in [31, 32], reduces the description of the long-time 
behaviour of (4.5) to a spectral problem for an auxiliary operator K\(uj) defined below, 

(A) t = =^ [ e itz ((£*Y>o) fi Rl/ 9, (K x (u) - zr'AQ^o) dz. (4.6) 

JR-i 

The integral in (4.6) is understood as an improper Riemann integral, 

f f(z)dz = lim f f(x-i)dx. (4.7) 

JR-i a,b^oo J _ a 

We present here a heuristic derivation of (4.6), see however [31, 32] for a rigorous proof. 
First we get rid of the factor e~ ltLx in (4.5) by the following trick [23, 31, 32]: one can 
add to L\ a term which does not change the dynamics but which is s.t. the resulting 
sum K\ satisfies K\£lpfi = 0. The operator K\ is of the form 

K X = L + XI, 

for some non-self-adjoint operator / related to V, see (B.3). We present the explicit 
construction of K\ in Appendices A and B. It gives that e ltLx Ae~ ltLx = e ltKx Ae~ ltKx , 
and that e -1 *^^^ = £|a,o- Using these relations we obtain from (4.5) 

(A) t = ((BVt>) ® ^R,/3, e^M^o} • (4.8) 

The operator K\ is unbounded and non-self-adjoint, and we do not know a priori 
whether e ltKx is defined on A£lp . Consequently, the formula above and the next one 
are formal, heuristic expressions (which obtain a rigorous meaning using a complex 
deformation, see below and [31, 32]). Furthermore, we represent the propagator as an 
integral over the resolvent, [31, 32] to arrive at 

(A) t = / e" 2 ((B*fo) n Rjf3 , (K\ - z^Anpp) dz. (4.9) 

lm JR-i 

Next we "uncover the resonances". We use the following notation for the second 
quantization of a a one-body operator O acting on single-particle wave functions of the 
variable k G M. 3 : 

dT(0)= f a*(k)Oa(k)d 3 k. (4.10) 
Jr 3 

Let 

N = dr(i) ® 1 R + 1 R ® dT(l) 

be the total number operator on Tin = ^"(L 2 (IR 3 , d 3 k)) T(L 2 (R 3 , d 3 k)), and define 
the operator 

D = dr(t?)®l-l R (8)dr(<?), (4.11) 

where ■& = ^(k ■ V + V ■ k), with k = . The operator D is self-adjoint (even though 
dr(i?) alone is only symmetric but not self-adjoint, see [31]). Let 

U(u) = e- iu}D , (4.12) 
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lo £ R, be the group of unitary transformations on Tin generated by D. We transform 
the operator K x unitarily as K x (uj) = U{uj)K x U(ujY 1 . From -i[dr(tf), dr(|&|)] = 
dr(l), and from equations (2.20), (2.21) and (2.27), we obtain the relation i[D, Lq] = N. 
Consequently, 

K x (u) = L + ujN + XI{uj), (4.13) 

where I(uj) = U (uj)IU (uj)" 1 . The explicit formula is presented in Appendices A and 
B. Since U(u) is unitary (for to G R), and since both Atlp^ and (B*ip ) (g> ^ are 
invariant under U(oj), the r.h.s. of equation (4.9) equals that of (4.6), for all uj £ R. 
The integrand in (4.6) can be extended to an analytic function in uj, for w in a strip 
< Imcj < ujq = ^ (which is continuous as Imuj J, 0). Since the r.h.s. of (4.6) is 
constant in uj for real to, it follows that (4.6) holds for all complex uj in the strip. This 
concludes our heuristic derivation of formula (4.6). 

Let us take to = ico', for some uj' > 0. The point is that the spectrum of K\(uo) is 
much easier to analyze than that of K\. This is so because the spectrum of Kq{lo) = 
L + \JN is 

spec(KoM) = {{Ei - Ej} iJ=ltmmmtN ) U n >i (icu'n + R). (4.14) 

The eigenvalues Ei — Ej have eigenvectors ipi (g) ipj ® ^r,/? and the lines io/n + R are 
horizontal branches of continuous spectrum. There is a gap of size to' separating the 
eigenvalues from the continuous spectrum of Kq(lo), so if we add the perturbation 
XI (u>), which is bounded relative to Kq(u>), we can follow the location of eigenvalues by 
simple (analytic) perturbation theory, provided A is small compared to oj'. One obtains 
the following result. 

Theorem 4.1 Fix lo' > 0. There is a constant cq > s.t. if |A| < cq/ (5 then, for all 
uj with \mu) > uj' , the spectrum of K\(uj) in the complex half-plane {Imz < uj'/2} is 
independent of uj and consists purely of the distinct eigenvalues 

{4 s) (A) | e€spec(L s ),s = l,...,i/(e)}, 
where 1 < v(e) < mult(e) counts the splitting of the eigenvalue e. Moreover, we have 

(s) (s) 

liniA^o ke (A) — e| =0 for all s = 1, . . . , v(e), and we have lme e (A) > 0. Also, the 
continuous spectrum of K\(uj) lies in the region {Imz > 3c<//4}. 

Remarks. 1. The proof of the theorem uses standard analytic perturbation theory. 
The constant Co in Theorem 4.1 is defined by the condition that the eigenvalues ei s \\) 
stay away from the continuous spectrum of K\(uj). Since 

(a) the continuous spectrum of K\(uj) lies in a neigbourhood of order A around the 
continuous spectrum of Kq(uj), i.e., in the region {Imz > 3u//4}, provided |A| < c±uj' , 
for some constant c\ (see (4.14)), and 
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(b) the eigenvalues e e (A) have imaginary part bounded from above by C2A 2 , for 
some constant C2, see (1.5), 

the eigenvalues and the continuous spectrum of K\{w) are separated provided |A| < 
min{ci,C2}w / < 27rmin{ci,C2}//3. 

The fact that the spectrum of K\{uj) must lie in the closed upper half plane is quite 
clear since (A) t , (4.9), must stay bounded as t — ► 00. 

2. By construction we have K\(u)£lpfi = 0, so we set = 0. 

We now perform the pole approximation by deforming the contour z = R — i in the 
integral (4.6) into the contour z = R + ^[/x + u//2], where we introduce 

/i = max{Im4 s) (A)}. (4.15) 

e,s 

This contour separates the eigenvalues (A) (lying below it) from the continuous spec- 
trum (lying above it). Indeed, the gap between the eigenvalue with biggest imaginary 
part and the continuous spectrum is bounded from below by 3u//4 — ji. 

With this contour deformation, we pick up the residues of the poles of the integrand, 
sitting at the resonance energies e^\\). Let be a small circle around not 
enclosing any other point of the spectrum of K\{u)). Then (4.6) gives 

(A) t = (4.16) 

-1 v(e) r 

^EE/ (S) e '* 2 <( 5 >o) ® Or./s, (K x (u) - zy'An^) dz + R, 

where 

R=^f j tz ((B*ip )®nR,(3,(Kx(u)-z)- 1 An /3 , )dz. (4.17) 

2TT1 Jr + ^ +u} '/2] 

We prove the following estimate on the remainder term R in Appendix D. 
Proposition 4.2 We have R = 0{\ 2 e'^ +u} ' ^). 

For simplicity of the exposition, we assume in the remainder of this section that 
the nonzero resonance energies / are simple poles of the resolvent {K\{u) — z)^ 1 
(i.e., that the ei^ are semisimple eigenvalues of K\(u); this will be satisfied in all our 
applications). It is easy to extend the following arguments to the general situation. 

Under the assumption above we can replace e by e e in the first term on the 
r.h.s. of (4.16), 

v(e) 

( A )t = E E ^ (( fl >o) ® n Rl/3 , Q^AQ^o) + 0(A 2 e-'*/ 2 ), (4.18) 

e s=l 
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where we introduced the (in general, non-orthogonal) projections 

= Q^(u) = ^ /" (#aM " *) _1 ds. 
If e{ s ^ is a simple eigenvalue of K\(w), then we have 

q£ s) = lx£ s) ><xS s) i, 

where the vectors Xe^ and Xe^ satisfy 

Kx(u)^ ] = ei s) xi s) and (A^))*xi s) = ^Fx« (4.19) 
and are normalized as 

(xi s \xi s) ) = l- (4-20) 

We obtain from (4.18) the relation 

((^))oo : = T lim \ f (A) t dt= ((B*i>o)®n R ^Q^Anw). (4.21) 

All the other terms vanish in the ergodic mean limit. 

If (A) t has a limit as t —> oo, as is the case when lmei s \\) > for ei s \\) / 0, 
then ((A)) ^ is just that limit (it may happen that (A) t does not have a limit, but 
((A)) oo always exists as we see from (4.21)). We have thus identified the limit term in 
the expansion (4.18) and we obtain the results of Theorem 3.1. 

For specific models, one can calculate (perturbatively in A, to any order) the reso- 

(s) (s) 

nance energies e e and the projection operators Q e , and one obtains estimates on the 
difference (A) t - ((A)) ^ 

Remark. As the results above show, the long time behaviour of averages (A) t 
is determined by the resonance energies ei^ and the resonance eigenvectors xi^ and 
Xe^ of the family K(ui) and its adjoint. This resonance data can be constructed from 
eigenvalues and eigenfunctions of operators A e a \ which act on f)s <8> f)s and which are 
independent of to. (This reconstruction can be viewed as a kind of "inversion" of the 
Feshbach map, see Sections 5 and 6.) We point out that this procedure is relatively 
easy to implement for translation-analytic systems, i.e., systems satisfying Condition 
A stated after (3.2). For systems wich do not possess this condition a renormalization 
group analysis in the spirit of [7, 8] has to be carried out, see also [31, 32]. 

5 Resonance energies ei\\) 

The goal of this section is to evaluate the main contribution (in A) to the resonance 
energies e*£\\), where e G {E m — E n | m, n = 1, . . . , A^}. Using a standard Feshbach- 
map argument (see e.g. [9, 5, 31, 32]) we obtain the expansion 

e W(A)= e -A 2 «W+0(A 4 ), (5.1) 
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(s) 

where the 5e £ C are the eigenvalues of the so-called level shift operator A e associated 
to e, defined by 

A e = P e IP e (L - e + iOj-^e/Pe, (5.2) 

where P e = P(Lg = e) <g> Pq is the eigen-projection of Lq associated to the eigenvalue 
e, P e = t — P e , and the restriction of P e L§P e to the range of P e is denoted by Lq (see 
also [29]). 

Consider an interaction of the form (2.32), where G is a symmetric matrix on 
f)S = C n . The following result follows from a direct calculation involving the explicit 
form of A e . We prove it in Appendix D. 

Proposition 5.1 (Level shift operators) The operator 

A e (e) := P e IP e (L - e + ie^PJPe 

has the representation 

A e (e) = A e , d (e) + A e , (e) + A e , m (e), (5.3) 

where the subscripts "d,o,m" stand for "diagonal, off-diagonal, mixed" and refer to the 
decomposition G = G d + G into a sum of a diagonal and an off- diagonal matrix in 
the energy basis. Let G be the matrix obtained by taking the complex conjugate of G, 
in the energy basis. We have 

2A e , d (e) (5.4) 
= P e {G A ®l-l®G d ){G d ®l+t®G d )P e (g,-J^g^ 

-p e (G d ®i-i®G d ) 2 p e (g jCO th(^ ~j^ g y 

2A e , (e) (5.5) 
= P e (G Q 1) [ U ]\ 9{H ^S (Ls -e + u + iey 1 (G Q ® t)P e 

+ F e (l®e-^G ) / U y {W \ U l (Ls-e + u + ie)' 1 (1 G Q e^)P e 

- Pe{G ®t)! U y ilul a JS (L s -e + u + i.)' 1 (1 e-f H *G a e$ B *)P e 

JRxS 2 I 1 - e H 

- P e (l®e-^-G e^) / ^M^j! ( Ls - e + u + ie )-i ( Go8l )P e 
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2A e , m (e) (5.6) 
= P e {(G d ® 1 + 11® G d ) , (G 1- 1® e-i Hs G e^ s )}p e x 

/ 

X \g' 2 i 2 

- P e {(G d ® 1- l®G d ) , (Go 11- 1® e-i Hs G e^ s )}p e 
x^,coth(^) ^^g), 
where {A, B} = AB + BA is the anti- commutator. 



Remarks. 1. If e is a simple eigenvalue then all factors e^2 Hs in (5.4)-(5.6) can 
be set equal to one. 

2. If the eigenvalues of H$ are non-degenerate then Ao(e) = Ao (e). In particular, 

(s) ' 

the eigenvalues Sq do not depend on the diagonal elements of G. 

To see that the statements in the remark are true, simply use in formulas (5.4)- 
(5.6), the facts that P = J2jPj ®Pj®PR, where pj = \ipj)(<Pj\ and P R = \tt Kj/3 )(n K>l3 \, 
and that P e = pi® pj <8> P R if e is simple. 

6 Resonance eigenvectors Xe^(A) and xi^(A) 

In this section we derive the following expansions of Xe^(A) and Xe^(A) (defined in 
(4.19)) in A: 



W 

e 



[1 - \P e {L Q {u) - e)- 1 P e /(^)P e ] # ® n R , p + 0(A 2 ), (6.1) 
[1 - XP e (L (u;) - e)- 1 P e (r)(o;)P e ] # ® fi R>/9 + 0(A 2 ), (6.2) 



(s) ~(s) 

where r/e <8> ^r,/3 and 7^ <8> ^ R /3 denote the eigenvectors of the level shift operator A e 
and its adjoint (A e )*, respectively, 

A e #®ft Ri/3 = # ® fi R>/ ,, (6.3) 
(A e )*#®ft R)/3 = ##®fi R/3 . (6.4) 

To arrive at the expansions (6.1), (6.2), we use the method of the Feshbach map (see 
e.g. [7, 8]), according to which we know that 

X W = (P e - \P e {K X {u) - S^r'PelH) Pe&\ (6.5) 

where K\{lo) = P e K\(uj)P e \ RanP e , and where G RanP e is the eigenvector of 

F Pe)£ ( a) (K A H) := P e (e - A 2 I(u;)P e (if A H - e^PJiu)) P e , (6.6) 
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with eigenvalue e e . We expand the r.h.s. of (6.6) in A using a Neumann series 
expansion of the resolvent P e (K x (u) - e { e s) )- l P e = P e (K {uj) + \I(uj) - ei s) (A)) -1 Pe- 
This gives 

F Pe £ ( s) (K A H) = P e (e - A 2 A e + 0(A 4 ))P e , 

where we also used (5.1). Consequently, taking into account the definition (6.3), we 
get the expansion 

#=#®^R,/3 + 0(A 2 ). 

Finally, a similar expansion of the r.h.s. of (6.5) gives (6.1). 

The analogous expansion (6.2) for xi^ is obtained by proceding as above, and by 
using the following little result. 

Lemma 6.1 The level shift operator associated to the eigenvalue e and the adjoint 
operator (K\(uj))* , Im.uj > 0, is given by 

P e (P)(uJ)P e (L (Lo) - e)- 1 P e (r)(^)P e = (A e )*. 

Proof. We have (recall (4.11) and (4.12)) 

(K x (u))* = [e- iujD K x e^ D ]* = e~' lWD Kl^ D = L (uJ) + A(/*)(cJ). 

Therefore we have the following expression for the level shift operator, 

P e (I*)(Zj)P e (L (uJ) - e^Pein^Pe 

= P e I*P e (L - e - iO)~ 1 P e rP e 

= [P e IPe(L - e + iOy^elPe]* 
= (Ae)*- 

In order to get rid of the parameter uj by analyticity, we have introduced in the first 
step the term — iO. ■ 

7 Special cases 

7.1 The qubit: proof of Theorem 3.4 

In this section we calculate the level shift operators and their spectral data for a qubit 
coupled to the Bose field. Recall that a qubit is a two-dimensional system, with state 
space (of pure states) f)s = C 2 , and Hamiltonian H$ = diag(£'i, E2). The interaction 
to the Bose field is given by (1.6) (which becomes (2.32) in the Hilbert space of positive 
temperature states). Recall the definition of the coupling parameter £(??), (1-8). We 
set 





' 1 " 




" " 


<P1 = 





, ^2 = 


1 



Let A = E2 — Ei > denote the gap in the energy spectrum of H^. The following 
result is an easy application of Proposition 5.1 to the specific model. 
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Proposition 7.1 In the basis {<pi ® (pi,f2 ® ¥2}, the level shift operator Aq is given 
by 

e -/3A/2 _ x 



9A ^2, ,2 £(A) 

zAn = — 17T C - 



(7.1) 



cosh(/?A/2) [ -1 e^ A / 2 
27ie (one- dimensional) level shift operator Aa is given fry 

2A A = (b 2 -a 2 ){g,u- 1 g) + \c\ 2 VN. f u 2 | 5 (|u|, a)| 2 coth(^) * 

-i7r 2 |c| 2 ^(A)-i7r(6-a) 2 e(0), (7.2) 



and A_a is obtained from the r.h.s. of (7.2) by switching the sign of the real part, 
A-a = -Aa. 

Remarks. 1. The vector ^s,/3 V?i ® Vi + e~@ A / 2 (f2 <8) 9?2 spans the kernel of Ao. The 
gap in the spectrum of Ao is 7r 2 |c| 2 £(A), it is exactly twice the gap of ImA-tA coming 
from a G with constant diagonal in the energy basis. 

2. £(0) is non-zero only if the infra-red behaviour of g is given by p = —1/2, c.f. 
(1.7). 

It is also easy to obtain the eigenvectors of the level shift operators (and their 
adjoints, see Lemma 6.1). 

Proposition 7.2 For the qubit model, the vectors determined by (6.3) and (6.4) are 

rii 1] = # } = = (1 + e-^)- 1 / 2 ® Vl + e-^V 2 ^) , 
„( 2 ) = = (i + e -^)-V2 ( e -/3A/ 2v9i g, ^ _ ^ g, ^ f 

r]A = rjA = (p 2 ®<pi, 

V-A = f?-A = V?l ® V?2, 



f/iey are normalized as (ry e | = 1, as required by (4-20). 

Propositions 7.1 and 7.2 together with Theorem 3.1 imply Theorem 3.4. ■ 

Let us show the assertion of Remark 5 after (1.11). Suppose that /5o = \ipj)((fj\ for 
a fixed j = 1,2. This initial state is represented on Tis = C 2 ® C 2 by the vector 
tpo = ipj®<Pj. The operator B defined by relation (4.3) is B = Cj~&s®\<Pj}(<Pj\, where Cj 
is a complex number. Thus B*i/jq = Cjipo and we may replace B* by Cj in formula (3.3), 
Theorem 3.1. By Proposition 7.2 we have Q±a = |?7±a)(?7±a| <8> |^R„a)(^R,/?| + 0(\ 2 ), 
so the terms in (3.3) coming from the nonzero resonances are 

cj (tp ® ^r,/3, Q±AAttfr ) = cj (ipj ® ipj ® n Rj/9 , Q±AAnp >0 ) = 0(A 2 ), 

since [<fj <X> <Pj,?7±a) = 0. The calculation of the constants Co for j = 1,2 is carried 
out in the same fashion. Here one uses that A = p\ t i ® Ig <8> 1r and the explicit form 
of ^,0, see (4.4) and (2.8). 
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Finally we show how to arrive at the relations given in the illustration involving 
(1.15). As one easily verifies the state Vo representing (1.15) in f)g <g> f)g is given by 
ipo = -^(fi <8> f2 + f2 <8> ¥2)- The associated operator B (see (4.3)) is 





e /3E 1 /2 e /3E 2 /2 

The relations after (1.15) follow by easy direct calculation, using Theorem 3.1 and 
Propositions 7.1 and 7.2. 

7.2 Energy conserving interactions 

In this section we apply the dynamical resonance theory to interactions which com- 
mute with Hs- Such quantum non-demolition interactions are widely studied in the 
literature, see e.g. [15, 17, 35, 33, 40] and references therein (this is a small fraction 
of works on this subject - there is an immense number of papers on it). Some energy 
conserving models can be solved explicitly and provide benchmark examples for our 
methods, as well as starting points for a perturbative treatment (for small [Hs,v]). 

Let f)s = C N , Hs = diag(£i, . . . ,E N ), where Ei < E 2 < ■ ■ ■ < E N . (A similar 
analysis can be carried out if there are degenerate eigenvalues.) We couple S to the 
bosonic reservoir via the interaction Xv, (1.6), which is represented on the positive 
temperature Hilbert space by the operator XV of the form (2.32), with a matrix G that 
commutes with Hs, i.e., 

G = diag(7i, . . . ,7^), 

with 7j £ K. Let d > 1 be the dimension of the quantum field, and uj = \k\. The 
dynamical resonance theory yields the following result (which we prove in Appendix 
D). 

Proposition 7.3 Theorem 3.1 implies that [~p t ] m ,m = \po]m,m, o,nd that, for m / n, 

[ptUn = ([Pokn + 0(A 2 )) e --AEm-E n)+ iXH [ 5 En . Em+ 0^) ] + Q{> 2 ^ /2^ (? _ 3) 
where 



( ifp > ^ 



5 En -E m = \{l 2 m - 7n) (9,^ 1 g) + Kim ~ In? < 



<(o)>o if v = if ( 7 - 4 ) 

+00 ifp < 



Here, p is the power characterizing the infrared behaviour of g, c.f. (1-7), and £(0) is 
given by (1.8) in which the integral is taken overM. d instead ofM 3 . 

As it turns out, this model is explicitly solvable. For illustration purposes we com- 
pare the explicit solution to the results obtained by our dynamical resonance method. 
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Proposition 7.4 (Explicit solution) The reduced density matrix elements are given 
by 



and where u(k) = \k\, a(k) € S d 1 . 

Remarks. 1) The same result for the decoherence function T(t) has been obtained 
in [35, 33]. There, the strategy was to obtain an explicit expression for the matrix 
elements, first for the system where the reservoir's modes are discrete, and to take 
afterwards the continuous mode limit. Our calculation is performed directly on the 
system with continuous modes. We give it in Appendix D. 

2) We have (g,oj~ 1 g^ < oo, as follows from the assumption fjbpg G L 2 (M. d , d d k). The 
latter is required for the interaction (1.6), (2.33) to be defined (recall that \ip has an 
^^-singularity at zero, see (2.17)). 

The next statement is evident from the explicit form of the matrix elements. 

Corollary 7.5 The populations are constant, \p t ]m,m = [Po]m,m f or a H m an d M t. 
Full decoherence takes place if and only ifT(t) — ► oo as t — > oo, i.e., if and only if 

The following result examines the asymptotic behaviour of a mj „ (t) . It illustrates the 
compatibility between the explicit solution and the result obtained from the resonance 
theory. Its proof follows readily from the explicit expressions given in Proposition 7.4. 

Proposition 7.6 Let 5E n -E m be defined as in (7.4). We have 



For p = the decoherence function grows asymptotically linearly in time, T(t) ~ 
7r^(0)t, for large t. 

Remark: Comparison with the results in [35]. Propositions 7.3 and 7.6 shows the 
following behaviour of the decoherence function. In dimension d = 1, for an infra-red 
behaviour of the coupling function g(k) ~ |/c| 1//2 , and for large times, the decoherence 
function is T(t) ~ £(0)i, which becomes T(t) ~ Tt, for small temperatures (see Remark 
11 after (1.11)). This is the same result as obtained in [35]. In dimension d = 3, 
the contribution to the decoherence function which is quadratic in A vanishes, see 
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Proposition 7.3. This is in accordance with the result of [35] that decoherence is 
incomplete in three dimensions. 

We conclude this section with the analysis of the resonances for energy-conserving 
interactions. The next result (see Appendix D for a proof) shows that the kernel of 
Lq is invariant under energy-conserving interactions. In particular, the degeneracy of 
the eigenvalue zero is not lifted under perturbation. This means that there are no 
resonances bifurcating out of the origin. 

Proposition 7.7 (Zero resonances) The kernels of Lq, K\ and K\(uj) coincide. 
They are spanned by the vectors 

is) 

Xq = ¥s <8> fs <8> ^R,/3, s = l,...,N. 

The kernel of (K\(lj))* is spanned by 

Xq ] = <Ps <8> ¥s ® (1 + T s )Q Rj p, s = l,...,N, 

where T s = 0(A) satisfies (Or,/?, T s tt R> p) = 0{\ 2 ). In particular, there are no reso- 
nances bifurcating from the origin. 

As an application of this result, we show, using resonance theory, that the popula- 
tions are constant (a fact we know already from Corollary 7.5). Proposition 7.7 gives 
us the following expresssion for the projection onto the zero resonances: 

Q = f (K x (uj) - z^dz (7.5) 

N 

= ^^Ifs ®<Ps® ^R,/3>(^s ® ® (1 + T s )Q Rjl3 \, 
s=l 

where 

Ws = ((i + T s )n R , p , n R , p ) = i + (T;) n ^ = 1 + o(a 2 ) 

is a normalization factor. Using (7.5) in equation (4.21) we obtain 

e -/3E m /2 

((Pn,m))oc = <Wn 7= {B*1p , (f m <g> <p m ) 

V Z S,/3 

= 0~n,m ((B*1p ) <8> ^R,/3, (p n ,rn ® % ® 1r)%o) 

= 0~n,m {Pn,m)f—Q • (^-6) 

To see how to recover the fact that the populations are constant in time we note that 
since p ntn ® Is <8> 1r commutes with K\(u>) (this is not true for p n>m , m / n) and since 
^(cjjfi^o = we have 

(K X (UJ) - z)- 1 (p ntn ® Is ® 1r)%0 = -2 _1 (Pn,n ® Is ® 1r)%0- 
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The contour integrals in (4.16) can then be evaluated, 

1 f e itz 

(Pn,n) t — (Pn,n) t=0 = (Pn,n) t=0 ^"7 / — — dz. 

The last integral has to be understood in the sense of the appropriate limit (4.7). By a 
standard contour deformation we see that we may take u/ > in this integral as large 
as we please, so the r.h.s. is zero (see the proof of Proposition 4.2 in Appendix D for 
details). Consequently, (p n ,n) t = (Pn,n) t=Q , i.e., [p t ]n,n = [Po\n,n- 

Next we examine the nonzero resonances. We have proven above that, under energy- 
conserving perturbations, there are no resonances bifurcating out of the origin. Now we 
show that nonzero eigenvalues may migrate into the upper complex plane or may stay 
on the real line, depening on the energy-conserving interaction. We know from (7.6) 
that the ergodic means of the off-diagonal matrix elements tend to zero as t — ► oo. The 
off-diagonal matrix elements can be purely oscillatory, which corresponds to resonances 
staying on the real axis, or they decay in case the resonances move into the upper 
complex plane. To examine which of the two cases happens we consider the level shift 
operator A e associated to an eigenvalue e of Lg, given in (5.3). Since G is diagonal, only 
the diagonal term (5.4) is present. Let I e be the set of indices m,n s.t. e = E m — E n , 
where the Ej are the eigenvalues of Hs- We have P e = J2i e Pm'8 l Pn, where pj = \ipj)(<Pj\, 
and we obtain easily the relation 

A e = J2 Pm ® Pn [(7m ~ Tn) (9,^9) ~ i<(0)(7m - Inf] , 
h 

see also (1.8) (where the integral is understood to be taken over M. d ). This shows that 
if the coupling to the low energy modes of the reservoir is weak, p > then the 
imaginary part of all the level shift operators vanish, since £(0) = 0. Note that for e = 
the sum extends over p m ®Pm and hence ImAo = 0, which is consistent with the result 
of Proposition 7.4. Furthermore, if G is constant on RanP e then ImA e = as well. In 
these cases the resonances stay on the real line, to second order in A. Proposition 7.4 
asserts that this actually happens to all orders in A, a result which one could extract 
from the resonance theory as well, by performing a deeper analysis. For instance, one 
sees from the explicit form of the interaction operator I(u>) (Appendix B, after (B.15)) 
that 

I(u)(p m ®lf n ® £Ir : {3 = (7m ~ 7n)(ls ® % <8> U(u)(pp(g)) <fi m ® H>n <8> ^K,f3, (7.7) 

and consequently, that if 7 m = j n , then the Feshbach map applied to K\(uj) is simply 
equal to eP e . The isospectrality of the Feshbach map implies that the eigenvalue e of 
Lq does not move under the perturbation (to any order in A). 

7.3 Observations about quantum registers 

Consider a register of L qubits located at positions xj, j = 1,...,L, which do not 
interact directly with each other. This L-bit register is placed in an environment, 
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modelled by a thermal Bose field described in Section 2. It is known [35, 17] that 
the decoherence properties of the register depend on the relation between the spacing 
of the qubits in the register and the correlation length of the reservoir. Consider the 
following two cases: 

(i) The qubits are very far apart, s.t. min^y |xj — Xj\ is much larger than the 
correlation length of the reservoir. 

(ii) The qubits are packed very closely together, so that much 
smaller than the correlation length of the reservoir. 

In case (i) we expect that a good approximation to the true dymamics is given by a 
register where each qubit is coupled to its own reservoir. Such a system is described 
by the zero-temperature Hilbert space of pure states 

L 
3=1 

where hs,j and f)Rj are the Hilbert spaces of the j'-th qubit and the j-th reservoir. 
(Recall the notation of Sections 1 and 2.) The subindex i indicates that the qubits are 
coupled to individual, independent reservoirs. The Hamiltonian is given by 

L 

H hX = J2Hx,j, (7.8) 

3=1 

where H\j acts nontrivially only on the j-th qubit-reservoir pair, and is given by (1.1). 
Accordingly, the positive temperature Hilbert space is Hi = ®j =1 Hs,j <8> Hrj, and 
the Liouville operator has the form L\ = J2j=i where each L\j generates the 
dynamics of a single qubit coupled to its own environment, as in Section 3. Let us 
label the state <p mi <8) • • ■ <8> i Pm h of the register by m = (mi, . . . , uil) € {0, 1} L . The 
reduced density matrix for the register, p t , is simply the product of the reduced density 
matrices of the single qubits, ~pj t . Its matrix elements are 

L 

[/?t]m,n = JJlPj'jtlm^nj) 

3=1 

where each [pj^m^nj evolves according to (1.9)-(1.11). 

Next let us consider the case (ii). We expect that the true dymamics is well ap- 
proximated by an interaction term where all the qubits sit in the same location. In this 
situation one observes collective decoherence, [35, 17]. Here, the subindex "c" stands 
for "collective" . The Hilbert space of pure states (at zero temperature) is 
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where f)sj is the Hilbert spaces of the j-th qubit and f)R is that of the reservoir. The 
Hamiltonian is given by 

H^^Hsj + Hn + xl^GA ®v(g), (7.9) 

3=1 \3=1 J 

where H^j is as above and where Gj acts as a fixed matrix G on the j-th qubit and 
trivially on all other qubits. The positive temperature Hilbert space is 



where Hj = C 2 <8> C 2 , and Hr is given in (A.l). The generator of dynamics takes the 
form L c \ = Lq + XV C , with 

L o = Y + Lr and Vc = f Y G * ® % J Ms) ( 7 - 10 ) 

(compare also with (3.1)-(2.33) and (2.29)). Here, Lgj acts non-trivially, as i^s®!— 
-ffS) only on the j-th qubit space, and Gj (8> lj acts as G<8> 1 on the j-th qubit space and 
trivially on the other qubits. Our general result, Theorems 4.1 and 3.1 and Proposition 
5.1 are valid for this model. However, the analysis of the level shift operators associated 
to the Liouville operator (7.10) becomes increasingly more involved with growing L. 
This is simply due to the size of the matrices representing the level shift operators. 
We point out, though, that the particular structure of V c may facilitate the spectral 
analysis of the level shift operators. Furthermore, since we have explicit formulas for 
the levels shift operators, our method may be suitable for a computer-based analysis. 

The Hamiltonians and H c \, (7.8) and (7.9), are in a certain sense extreme 
cases, as described by (i) and (ii) above, of an intermediate Hamiltonian. The latter is 
given in [35] as 

L L 

H x = H s ,j ® 1 + 1 ® #r + A Gj <g> v? (e-**^) . (7.11) 

3=1 3=1 

Here, g is a fixed form factor, typically imagined to be so that its inverse Fourier 
transform, g, is peaked around the origin in x space. The Fourier transform of the 
function g shifted to the position of the j-th qubit, gj(x) := g(x — Xj), is then just the 
e~ lkx ig appearing in the interaction of (7.11). 

The following are two problems in connection with quantum registers: 

PI Derive the form of Hamiltonian (7.11) from a real physical situation, e.g. where 
qubit j is represented by two levels of an atom located at position Xj. Quantify 
the approximation schemes regarding cases (i), (ii) above. 
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P2 Analyze the level shift operators associated with the interaction (7.10), analyti- 
cally or numerically. Obtain the dynamics of reduced density matrix elements (as 
in (1.9)-(1.11)) for the collectively decohering quantum register. Identify coherent 
subspaces. 

A Unitary transformation of the positive-temperature 
Hilbert space 

It is convenient to work with a unitarily transformed version of the Araki- Woods rep- 
resentation for the Bose field (2.14)-(2.18), c.f. [22, 18]. The Hilbert space (2.14) is 
transformed unitarily as 

F(L 2 (R 3 , d 3 A;)) ® F{L 2 (R 3 , d 3 k)) ^T = F(L 2 (R x S 2 , duda)), (A.l) 
according to the map 

a*(f 1 )---a*(f m )n®a*(g 1 )---a*(g n )n 

— a*(x+T/i) • • • a*{ X+ Tf m )a*(x-Tgi) ■ ■ ■ a*{ X -Tg n )Sl, (A.2) 

where the vectors Vt are the vacua in the respective Fock spaces and, and where the 
a# are the creation and annihilation operators in the respective Fock spaces. We have 
introduced T which maps functions f(k) 6 L 2 (R 3 ,d 3 k) into functions (Tf)(u,a) € 
L 2 (R x S 2 ,duda), according to 

\ f f( u i a ) if U > 0, , . „. 

(Tf)(u,a) = u < is-r, \ ., n ( A - 3 ) 
v J/K ' ' \ -e 10 /(-n, a) if u < 0, y ' 

where / is represented in polar coordinates and <p is an arbitrary real phase. This 
phase is a parameter which can be chosen appropriately to satisfy Condition (A) after 
(3.2) for a given coupling function g, see [18]. The x± m (A.2) are indicator functions, 
X+(u) = 1 if u > 0, x+(u) = if u < 0, and x- = 1 - X+- 

One verifies that the thermal annihilation operators, represented in the Araki- 
Woods representation by (2.15), take the following form in the unitarily transformed 
system: 

ap{f) = a{\J 1 + »p(u)x+(u)uf(u, a)) - a* ( e 1 ^ n p {-u)x- {u)uf(-u, a)) . (A.4) 

(The a*p(f) are obtained by taking the adjoint on the r.h.s. of (A.4)). A short calcula- 
tion shows that the thermal field operator (2.33) becomes, via the unitary transforma- 
tion, 

for / G L 2 (1R 3 ), where fp is given in (3.2), and where the ip in the r.h.s. is the field 
operator in T . The equilibrium state is represented by the vacuum vector of J 7 , 

n R ,f3 = Q. (A.6) 
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For a one-body operator O acting on wave functions of the variables (u, a), we write 

dT(O) = ! a*(u,a)Oa(u,a)duda. (A.7) 

for the second quantization of the operator O. The new representation has the advan- 
tage that the dynamics of the field is generated simply by 

L R = dl», (A.8) 

the second quantization of the operator of multiplication by u. We have LrOr,^ = 0, 
and for z G C, 

e^V/3(/)e-^ R = 2" 1 / 2 Me"* 4 /) +^(e zu f)) , (A.9) 

which gives the dynamics for z = it. 

It follows from (2.29), (2.32) and (A. 5) that the Liouville operator L\ acting on 
Hs ® T is given by 

L\ = L + XV, (A. 10) 

L = L s + L R = i^s®lls-ls^^s + dr(u), (A.ll) 
V = G® i s ®<p(g(i). (A.12) 

One verifies that the operator D, (4.11), is represented in the unitarily transformed 
space as dT(id u ), it generates translations in the variable u G R (see also [31]). The 
unitary group U(u), (4.12), is thus given by the translation group 

U(uj) = e ~^r(id u )_ (A. 13) 

The spectrally deformed Liouville operator acting on Hr, (g> T is 

L A (w) = L + L,jiV + Ay(w), (A.14) 

where = dr(l) is the number operator in J 7 , and where V{uj) = e -" dr ( 9 u)y e ^dT(d u ) 
(see also (B.14)). 

Observables of R are operators that can be built from sums and products of thermal 
creation and annihilation operators. Strictly speaking, one considers bounded operators 
built from exponentiated thermal field operators, e 1 ^^. They form the so-called 
Weyl-algebra, a von Neumann algebra 

m R C B{F), (A. 15) 

see e.g. [13]. 



36 



B The operators K\ and K\{u) 



We consider the positive temperature Hilbert space in its form given in Appendix A. 
The operator K\ can be expressed in terms of the non- interacting Liouville operator Lo, 
the interaction V, see (A.10)-(A.12), and modular data J, A (see e.g. [13]) associated to 
the vector £lp t o and the von Neumann algebra 971, (A. 15). J is an anti-unitary operator 
and A is a self-adjoint non-negative operator. The defining properties of J and A are 

JA^Mfiflo = M*Q Pfi , (B.l) 

for any M G 97T, where M* is the adjoint operator of M. From this property and the 
facts that Lo^/3,o = and V = V* we readily see that the operator 

K x = L + XI, (B.2) 
I = V - JA 1/2 VJA 1/2 , (B.3) 

satisfies K\Qqq = 0- The operators J and A satisfy JA = A _1 J, and therefore 
JA 1 / 2 VJA 1 / 2 = JA 1 I 2 VA~ 1 / 2 J. The theory of von Neumann algebras (Tomita- 
Takesaki) tells us that conjugation with A 1 / 2 leaves the algebra 971 invariant (provided 
the operators in question exist; A is unbounded), and furthermore, that all operators 
of the form JMJ, M G 971, commute with all operators G 971, so the subtracted term 
in (B.2) commutes with all observables of S + R and hence does not alter the dynamics. 

We now give explicit expressions for the operators J, A and K\, see also [13, 18, 
31, 32]. The modular data is 

J = J s <g> Jr and A = A s <8> A R , (B.4) 

where 

A s = e _/3Ls , (B.5) 

Ar = e-^, (B.6) 

J S (f>l <S) (j)r = C(p r ®C(j) U (B.7) 

■/R^nO^c-i, ■ ■ ■ ,u n ,a n ) = e m ^ n (-ui,cJi,...,-ti n ,cr n ), (B.8) 

where the action of the antilinear operator C is to take the complex conjugate of vector 
coordinates in the basis {ipj}^ =l of Hs, and ip n is the complex conjugate of ip n € T 
(see (A.l)). The phase <f> G M is the one appearing in (3.2). It may be chosen suitably 
to satisfy condition (A) (see after (3.2)), given a form factor g. Relation (B.8) shows 
that 

J R a*(f(u,a))J R = a*(e^J(-u,a)), (B.9) 

for / G L 2 (R x S 2 ). 

We use relations (2.32) and (B.4)-(B.7) to obtain 

I = V-V, (B.10) 
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where 

V' = l s ® e~% Hs Ge% Hs ® 

v2 



a*(^)+a(e-% /3 )l. (B.ll) 



We have set G = CGC here, and we recall that gp is defined in (3.2). 

We now give the explicit form of the spectrally deformed operator K\{u) = Lq + 
luN + I(uj), where U(u) = e ~ lujdr ( l9 ^) (see also (A. 14)). The transformation under 
U (u) of creation and annihilation operators is given by 

U(io)a*{f)U(ijy 1 = a*{f{- + u)), w£R, (B.12) 

where /(• + oo) is the shifted function (u,a) f(u + u>,a). Relation (B.12) can be 
written in the form U{ui)a^ (f)U(uj)^ 1 = <x&(& )du /). In order to obtain an analytic 
extension of (B.12) to complex u, we need to take the complex conjugate of co in the 
argument of the annihilation operator (since the latter is anti- linear in its argument). 
We thus have (see also (A. 5)) 

I(u) = V(u)-V'{u), (B.13) 
V(uj) = G®ls®-^=K( 5/3 (-+ W )) + a(^(- + cJ))], (B.14) 

V'{w) = Is® e-^ s Ge^®I \a*(gp(- +a(e-M^gp(- . (B.15) 

V2 L -I 

Let us finally prove the validity of (7.7). Definition (A. 13) of U(u>) implies that 

Further, it follows from (B.5)-(B.8) that 

Therefore, we have 

V'(oj)<p m ® V9„ (8) ft Ri /3 

= e-^ Em - En )l 2 U{u) Jk 1/2 V^ n ®^ m fi R>/9 

= 7 ™e- /3(i?m - £ " )/2 ?7(^)JA 1 / 2 (lls ® Is ® M<?)K 8^ m ® ^ R/3 

= 7 n [/(u;)(]l s l s tpp(g))<p m <8> <£> n <8> R)/3 , (B.16) 

where we have also used in the third step that Gip n = jn^Pn, an d in the last step that 
J R A R <^g(<7)J7 Rj/ g = <pp(g)QR t i3 (which follows from (B.l) and the fact that <pp(g) is 
self-adjoint). Combining (B.16) with 

V(uj)cp m (g> ip n <g> fi Rj; g = 7 m f7(u;)(ls <8> Is ® tp/3(g))tPm ® ¥>n^R,/3 

and with the fact that J(u;) = V(u) — V'(u), we obtain (7.7). 
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C Perturbation of equilibrium states and outline of proof 
of Theorem 3.3 



We give an outline of the expansion of the equilibrium state p(/3, A) of a coupled system 
(or its reduction to a subsystem) in terms of the coupling constant A. Explicit calcu- 
lations for the situation of the qubit yield a proof of Theorem 3.3. We do not carry 
them out here. 

Consider an TV-level system coupled to the Bosc field, as described in Sections 2, 3. 
It is well known [13, 9, 18, 22, 29, 31, 32] that the equilibrium state A) w.r.t. the 
interacting dynamics, (2.28), is represented by the vector 

Here, the operator L\ and the equilibrium state of the uncoupled system, fi/^o, are 
given by (2.29) and (4.4), respectively, and TL is the Hilbert space (2.25). Our task is 
to expand the average 

[Poo\m,n = <ty?,A, (Pn,m ® ls)tyj,A> (C2) 

in powers of A, where we recall the definition p n ,m — Wn)( i Pm\-, the {ipj} forming an 
orthonormal basis diagonalizing the Hamiltonian Hg. Since Lo^/3,o = 0> it is convenient 
to use the Dyson series, 

e-^/ 2 % = e-^/V L °/ 2 %, 

rP/2 r s n _i 



J±L fP/i rs„-i 

fyj.o + T"(-A) n / dai ■ ■ ■ / ds n V( Sl ) ■ ■ ■ V(s n )n p>0 , 

n=l Jo J ° 



where V(s) := e sL oy e sL °. Accordingly, it is clear how to arrive at an expansion of 
^/3,A ; (CI), and hence of the averages (C.2), in powers of A. ■ 

D Proofs of propositions 

Proof of Proposition 4.2. We expand the resolvent in (4.17) in A, 

(K x (oj) - z)- 1 = (K (u) - z)- 1 - X(K (u) - z)- 1 7(o;)(K A (a;) - z)~\ (D.l) 
The contribution to R coming from the free resolvent (Kq(lj) — z)^ 1 is given by 

^ [ . f(z)dz, (D.2) 

where f(z) := e ltz (B*i/}o, (Ls — z)~ 1 A^ls,f3), the inner product being that of 7^s. To 
arrive at (D.2) we use that A is an observable of S, that fi^o = £ls,p ® ^R,/3 and that 
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(Lr + uN)£lji t /3 = 0. The integral is understood in the sense of (4.7). Our first goal 
is to show that (D.2) is actually zero. Consider the integral J_ a f(x + + to' /2])dx. 
Since the only singularities of the integrand are poles on the real axis we can deform 
the contour of integration, yielding that for any r > ^[/j, + u//2], 

,-b 

f(x + i[fi + u>'/2])dx= (D.3) 

/b rr rr 

f(x + \r)dx+ I f(-a + iy)dy- f(b + \y)dy. 

-a J±\p,+u>'/2] ^i[/i+w'/2] 



/ 



It is easy to see that |/(— a + \y)\ < C e ty /a for a sufficiently large, and \f(b + iy)\ < 
Ce~ ty /b for b sufficiently large. It follows from (D.3) that 

fb 

f{x + ^ + u'/2})dx = f b _ a f(x + ir)dx + 0(l/a) + 0(1/6), (D.4) 



/ 



for a,b — > oo, and where the remainder terms are uniform in r. In the limit r — > oo 
the first integral on the r.h.s. of (D.4) vanishes, so by taking first r — ► oo and then 
a, b — > oo we see that (D.4) and thus (D.2) are both zero. 

Consequently, the term of order A in (D.l) does not contribute to the integral in 
(4.17). By iterating the resolvent equation we see that all terms with odd powers in A 
do not contribute either, because the interaction I is linear in creation and annihilation 
operators, and we take a "vacuum" expectation in (4.17). Thus we have 

R=^-r I e ii2 ((i?>o) ® fiR,/3, (K (u) - z)- 1 ^) x (D.5) 

x(K x (cu) - z)- 1 7(o;)(KoH - z)~ x An Pfi )dz. 

Finally we want to show that the last integral is 0{e~^ +u) Z 2 '). It is not hard to see 
that the norm of the integrand is bounded above by 

Ce-^ + "'M\\[(K Q (u;) - ^T^ill IK^oM - z)-^ 2 \\, 

for some constant C and where <J>i = (B*ipo) <g) JIr^ and <3?2 = AVLp^. In order to 
conclude that R = <3(A 2 e ~i^ +a/ / 2 ]) it thus suffices to show that 

f \\(K (u) - z)- 1 ^ j fdz < oo. (D.6) 

JR+^[H+ui'/2] 

The integrand in (D.6) is [(L - x) 2 + (u'N + to' /2]) 2 ]~ 1 ^ j ), which is read- 

ily seen to be integrable w.r.t. x € R (using for instance the spectral theorem 
for the commuting self-adjoint operators Lq and N). We have thus shown that 
= 0(A 2 e-i^'/2]). 

We point out that we did not use the specific form of 3>i in these estimates. The 
present argument works for all vectors <I>i corresponding to any initial state po on 
n s ® Hr. ■ 
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Proof of Proposition 5.1. This is an easy calculation using the following explicit 
form of all operators involved, as presented in Appendix A. 

Proof of Proposition 7.3. We refer to the paragraph after Proposition 7.7 for 
a proof of the fact that the populations are independent of time. We now concentrate 
on the off-diagonals. By using that G and H§ commute, it is not hard to see that the 
Feshbach map applied to K\(u), 

F Pe , z (K x (u)) = P e {e - \ 2 I(uj)P e {K x {uj) - z)- 1 P e /(w))P e , 

is diagonal in the basis {ip m ® ip n ® ^r,^} of RanP e , where m and n are indices s.t. 
e = E m — E n . By the reconstruction formula for eigenvectors, (6.5), it follows that 
the eigenvectors of K\(u) are of the form (p m ® ip n <8> (1 + r m n)fi^, where the T m ^ n 
are operators on Tin satisfying T m>n = O(A) and (Or,/?, T m ^ n Q. R ^) = 0(A 2 ). (Compare 
also with Proposition 7.7 and its proof.) In a similar way we see that the eigenvectors 
of [K\{uj)}* are of the form ip m ® tp n <& (1 + T m)fl )fiR )( g, for some operators T m ^ n on TL R 
satisfying f m , n = 0(A) and (sIrj, f m ^ K ^\ = 0(\ 2 ). Let d(e{ s) ) be the degeneracy 



of the resonance energy e[ s ^ . Then we have 

d(ei 3) ) 

Q ( e } = E <„,n s , (D.7) 

i=i 

(g) v?„ . 0(1 + 
where the normalization weights are 

Wm.jv.j ■= ((K + T mStj ,n s jn R ,p,(l + T mStj!ns jn R ^ = 1+0(A 2 ). 

We can now use expression (D.7) in expansion (3.3) to arrive at 

. , (s(m,n)) , , 

[PtUn = e lt£ ^-^C m>n ^o) + 0(X 2 e-^^ 2 ), (D.8) 

where 

= ^-^-JA 2 ^-^)^,^) (D.9) 
+iiA 2 ( 7 n - 7m ) 2 Hm ^, coth ^^g) + 0(A 4 ). 

The C m , n (ipo) in (D.8) are given by the following expression (see also the calculation 
leading to (7.6)) 



C m , n (^o) 

= {{B*ipo) <8> ^r,/?, y?„ (g> 9? m (g> (11 + r nim )Q Rj/3 ) 



e -/3E m /2 I , /^* \ 



y/~Zsj3 W n,m 

= [p ] m ,n(l + O(A 2 )). (D.10) 
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Relations (D.8), (D.9) and (D.10) show assertion (7.3) of the proposition, where 

$E n -E m ■= -^ 2 (ln ~ 7m) (a^^Q) 



+ 5 iA (7„ - j m ) lim ( g, coth 



Proof of Proposition 7.4. We absorb the coupling constant A into the matrix 
G (rescale G so that A = 1). The reduced density matrix is given by 

p t = Tr R [e- itL (p \n R ^){n Rt(3 \)e itL ] , 

where the trace is taken over Hr, (2.14) with L 2 (R 3 ,d 3 fc) replaced by L 2 (R d ,d d k), and 
where L = H$ + L R + G ® (p@{g) acts on the Hilbert space H = C N <8> Hr. Thus we 
have 

[PtUn = [PoUne-^-^ u^p (V^M^k) , (D.ll) 

where we denote the equilibrium state of R by lor^(-) = T^rH^r^) {^r,/3\ • )• Now we 
apply the Trotter product formula, 

UR,Je~'<™^ +L ^ e<™^ +L *A = lim lor, p ([X n (M)] M [X m (M)*} M ) , 

where X n (M) := e ~ itL ^ M e-^W^. Using that 

[X n (M)] M [X m (M)*] M 

= [X n {M)] M ~ l W p ^-A( 7n _ 7m)e --*/M^ [Xm (M)*] M -\ 

where Wp{f) := e'^^ is the thermal Weyl operator, and using the relation 

Wp{f)Wp{g) = e-^^Wpif + g), 
we obtain by induction the formula 

[X n {M)] M [X m {MT] M = exp j-i ^ 5^} W( 5M ), (D.12) 



where g M = -jj(ln- 7m) Ef=i e lfca; *f and 

,2 K-l 



t It 
S * = TH^n- 7m) (7n + 7m) ^ Im 



M 2 

fc=i 
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It follows from (D.12) and UR t p(Wp(f)) = e " 4</> coth (^/ 2 )/> ( see e .g. [13]) that in the 
limit M — > oo (where the sums over M turn into easy integrals which can be evaluated 
explicitly), we get 

lim ui K ,p{[X n (M)] M [X m (My} M ) 

= exp [-i( 7n - 7m )( 7n + lm )S{t) - ( 7n - 7m ) 2 r(t)] , (D.13) 

with Sit) and T(t) defined in Proposition 7.4. The proof of this proposition is now 
completed by combining (D.ll) and (D.13). ■ 

Proof of Proposition 7.7. We first notice that the Feshbach map (6.6) for e = 
and with spectral parameter z = vanishes, 

F Pofi (K x (u)) = -X^oI^PoiKxi^^PoI^Po = 0. (D.14) 

This is a simple consequence of the facts that RanPo is spanned by {ip s <8> ip s <8> Hr^} 
and that 

J(w) if s (g> 99 s ® fi Ri/ g = 0. (D.15) 

To see (D.15), simply use (B.2), (B.3), (B.10) and (B.ll), and that J s Ag /2 (G 
i s ) JsAg /2 = l s <8> G, which holds since [H s , G] = = [A s , G] and G is self-adjoint. 

Relation (D.14) implies that dimKeri^A = N. To obtain a basis for KevK\ we use 
the reconstruction formula (6.5) and notice that, again due to (D.15), only the term 
with Pq survives. 

Proceeding as in the proof of Lemma 6.1 it is readily seen that 

F P<hO ([K x (uj)}*) = [Fp <h0 (K x (u))}*=0, 
so dimKerf-fTA^)]* = N and we reconstruct a basis for Ker[Kx{oo)]* using (6.5), 

X ( S) = [P ~ APo([^A(w)]*) _1 Po[/(w)]*Po] ® fs ® ^R,/3. 

Since all the operators commute with the spectral projections of Ls we may "pull" 
ifs (g) if s through the operator [• • ■ ] to the left, and it is easy to identify the T s having 
the properties given in the proposition. Note though that equation (D.15) is not correct 
if we replace I{uj) by [I{uj)}*, so T s ^ 0. ■ 
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